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An approach is proposed to calculate Generalized Parton Distributions 



(GPDs) in a Constituent Quark Model (CQM) scenario, considering the con- 
stituent quarks as complex systems. The GPDs are obtained from the wave 



^ ■ functions of the non relativistic CQM of Isgur and Karl, convoluted with the 



GPDs of the constituent quarks themselves. The latter are modelled by using 
the structure functions of the constituent quark, the double distribution rep- 
resentation of GPDs, and a recently proposed phenomenological constituent 
quark form factor. The present approach permits to access a kinematical 
range corresponding to both the DGLAP and the ERBL regions, for small 
values of the momentum transfer and of the skewedness parameter. In this 
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kinematical region, the cross sections relevant to deeply virtual Compton scat- 
tering could be estimated by using the obtained GPDs. As an example, the 
leading twist, unpolarized GPD H has been calculated. Its general relations 
with the non relativistic definition of the electric form factor and with the 
leading twist unpolarized quark density are consistently recovered from our 
expressions. Further natural applications of the proposed approach are ad- 
dressed. 
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I. INTRODUCTION 



Generalized Parton Distributions (GPDs) [1] parametrize the non-perturbative hadron 
structure in hard exclusive processes (for comprehensive reviews, see, e.g., [2-7]). The 
measurement of GPDs would provide information which is usually encoded in both the 
elastic form factors and the usual Parton Distribution Functions (PDFs) and, at the same 
time, it would represent a unique way to access several crucial features of the structure 
of the nucleon [8,9]. By measuring GPDs, a test of the Angular Momentum Sum Rule of 
the proton [10] could be achieved for the first time, determining the quark orbital angular 
momentum contribution to the proton spin [9,11]. 

Besides, the possibility of obtaining, by means of GPDs measurements, information on 
the structure of the proton in the impact parameter [12,13] and position [14-16] spaces, is 
being presently discussed. 

Therefore, relevant experimental efforts to measure GPDs, by means of exclusive electron 
Deep Inelastic Scattering (DIS) off the proton, are likely to take place in the next few years 
[17-19]. 

In this scenario, it becomes urgent to produce theoretical predictions for the behavior 
of these quantities. Several calculations have been already performed by using different 
descriptions of hadron structure: bag models [20,21], soliton models [4,22], light-front [23] 
and Bethe Salpeter approaches [24], phenomenological estimates based on parametrizations 
of PDFs [25,26]. Besides, an impressive effort has been devoted to study the perturbative 
QCD evolution [27,28] of GPDs, and the GPDs at twist three accuracy [29]. 

Recently, calculations have been performed also in Constituent Quark Models (CQM) 
[30,31]. The CQM has a long story of successful predictions in low energy studies of the 
electromagnetic structure of the nucleon. In the high energy sector, in order to compare 
model predictions with data taken in DIS experiments, one has to evolve, according to 
perturbative QCD, the leading twist component of the physical structure functions obtained 
at the low momentum scale associated with the model, the so called "hadronic scale" , ji^. 
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Such a procedure, already addressed in [32,33], has proven successful in describing the 
gross features of standard PDFs by using different CQM (see, e.g., [34]). Similar expectations 
motivated the study of GPDs in Ref. [30]. In that paper, a simple formalism has been 
proposed to calculate the quark contribution to GPDs from any non relativistic or relativized 
model and, as an illustration, results obtained in the Isgur and Karl model [35] have been 
evolved from //q up to DIS scales, to NLO accuracy. In Ref. [31] the same quark contribution 
to GPDs has been evaluated, at /Xq, using the overlap representation of GPDs [6] in light- 
front dynamics, along the lines developed in [36] . 

In here, the procedure of Ref. [30] is extended and generalized. As a matter of fact, 
the approach of Ref. [30], when applied in the standard forward case, has been proven to 
reproduce the gross features of PDFs [34] but, in order to achieve a better agreement with 
data, it has to be improved. In a series of papers, it has been shown that unpolarized [37] 
and polarized [38] DIS data are consistent with a low energy scenario, dominated by com- 
plex constituent quarks inside the nucleon, defined through a scheme suggested by Altarelli, 
Cabibbo, Maiani and Petronzio (ACMP) [39], updated with modern phenomenological in- 
formation. The same idea has been recently applied to demonstrate the evidence of complex 
objects inside the nucleon [40], analyzing intermediate energy data of electron scattering off 
the proton. Besides, a similar scenario has been extensively used by other groups, start- 
ing form the concept of "valon", introduced more than twenty years ago [41] (for recent 
developments, see [42]). 

We here generalize our description of the forward case [37] to the calculation scheme of 
Ref. [30], in order to obtain more realistic predictions for the GPDs and, at the same time, 
explore kinematical regions not accessible before. 

In particular, the evaluation of the sea quark contribution becomes possible, so that 
GPDs could be calculated, in principle, in their full range of definition. Such an achieve- 
ment would permit to estimate the cross-sections which are relevant for actual GPDs mea- 
surements, providing us with an important tool for planning future experiments. Actually, 
as it will be shown, the proposed approach will be applied here in a Non Relativistic (NR) 
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framework, which allows one to evaluate the GPDs only for small values of the 4-momentum 
transfer, A 2 (corresponding to A 2 <C m 2 , where m is the constituent quark mass) and small 
values also for the skewedness parameter, £. The full kinematical range of definition of GPDs 
will be studied in a followup, introducing relativity in the scheme. 

The paper is structured as follows. After the definition of the main quantities of interest, 
an Impulse Approximation (IA) convolution formula for the current quark GPDs in terms 
of the constituent quark off-diagonal momentum distributions and constituent quark GPDs 
is derived in the third section. Then, the constituent quark GPDs are built in the fourth 
section, according to ACMP philosophy and using the Double Distribution (DD's) repre- 
sentation [3,26,43] of the GPDs. In the fifth section, as an illustration, results obtained by 
using CQM wave functions of the IK model and the obtained constituent quark GPDs will 
be shown. Conclusions will be drawn in the last section. 

II. QUARK MODEL CALCULATIONS OF GPDS 

We are interested in hard exclusive processes. The absorption of a high-energy virtual 
photon by a quark in a hadron target is followed by the emission of a particle to be later 
detected; finally, the interacting quark is reabsorbed back into the recoiling hadron. If the 
emitted and detected particle is, for example, a real photon, the so called Deeply Virtual 
Compton Scattering process takes place [8,9,11]. We adopt here the formalism used in Ref. 
[2]. Let us think to a nucleon target, with initial (final) momentum and helicity P(P') 
and s(s'), respectively. The GPDs H q (x, £, A 2 ) and E q (x, £, A 2 ) are defined through the 
expression 

1 1 irr^n A 

= H q (x, g, A*)-U(P>, s')iU(P, s) + E q (x, £ A 2 )-U(P>, s>) ^ » U(P, s) , (1) 

where A = P' — P is the 4-momentum transfer to the nucleon, ip q is the quark field and M 
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is the nucleoli mass. It is convenient to work in a system of coordinates where the photon 
4-momentum, = (qo,q), and P = (P + P')/2 are collinear along z. The £ variable in 
the arguments of the GPDs is the so called "skewedness" , parametrizing the asymmetry 
of the process. It is defined by the relation £ = — n - A/2, where n is a light-like 4- vector 
satisfying the condition n • P — 1. As explained in [9,11], GPDs describe the amplitude for 
finding a quark with momentum fraction x + £ (in the Infinite Momentum Frame) in a 
nucleon with momentum (1 and replacing it back into the nucleon with a momentum 

transfer A. Besides, when the quark longitudinal momentum fraction x of the average 
nucleon momentum P is less than — £, GPDs describe antiquarks; when it is larger than £, 
they describe quarks; when it is between — £ and £, they describe qq pairs. The first and 
second case are commonly referred to as DGLAP region and the third as ERBL region [2], 
following the pattern of evolution in the factorization scale. One should keep in mind that, 
besides the variables x, £ and A 2 explicitly shown, GPDs depend, as the standard PDFs, on 
the momentum scale Q 2 at which they are measured or calculated. For an easy presentation, 
this latter dependence will be omitted in the rest of the paper, unless specifically needed. 
The values of £ which are possible for a given value of A 2 are: 

< £ < v /z A 2 /v / 4M 2_r A 2 . (2) 

The well known natural constraints of H q (x, £, A 2 ) are: 

i) the so called "forward" or "diagonal" limit, P' = P, i.e., A 2 = £ = 0, where one 
recovers the usual PDFs 

H q (x,0,0)=q(x) ; (3) 

ii) the integration over x, yielding the contribution of the quark of flavor q to the Dirac 
form factor (f.f.) of the target: 

j dxH q { X) ^tf)=F?{tf)- (4) 

iii) the polynomiality property [2], involving higher moments of GPDs, according to 
which the x- integrals of x n H q and of x n E q are polynomials in £ of order n + 1. 
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In [30], the IA expression for H q (x, £, A 2 ), suitable to perform CQM calculations, has 
been obtained. 

Now, it will be shown that the same basic formula can be derived as the NR reduction of 
the definition (1) of GPDs, analyzed initially in the non-covariant framework of light-cone 
quantization, involving partons on their mass shell. 

Using the light-cone spinor definitions as given in the appendix B of [5], and defining: 
k + = (k + ks)/V2,k± = (ki,k 2 ), for the light-cone helicity combination s's = \\ = ++ one 
obtains 

Fl + (x, Z, A 2 ) = y/l^?H q (x, £, A 2 ) - -JL= Eq (x, £, A 2 ) , (5) 
so that, for £ 2 < 1: 

Fl + (x, £, A 2 ) = H q (x, Z, A 2 ) - e [\h q (x, Z, A 2 ) + E q {x, Z, A 2 )) + 0(£ 4 ) , (6) 

i.e. 

Fl + (x,Z,A*)=H q (x,Z,A 2 ) + 0(Z 2 )- (7) 

The reader should be aware that the so called "Munich Symmetry" for double distribu- 
tions excludes 0(Z) contributions to GPDs [44], so that the accuracy of the above equation 
is worse than it reads. 

According to the latter equation, in order to obtain the GPD H q (x, Z, A 2 ) for £ 2 <^ 1 one 
has to evaluate F\ +) starting from its definition, Eq (1). In the l.h.s. of the latter, using 
light-cone quantized quark fields, whose creation and annihilation operators, b\k) and b(k), 
obey the commutation relation 

{b(k'), tf(k)} = (27r) 3 2k + 5(k' + - k + )5 2 (k ± - k' ± ) , (8) 

and using properly normalized light-cone states 

(P'\P) = (2n f2P + 5(P' + - P + )5 2 (P ± - P' ± ) , (9) 

one obtains, for x > Z [2] 
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q.fet.A')-^/ ^ < t t( t + A) W > , (10) 



2P+V7 2^/|x 2 -e 2 |(27r)3 
where 

(tf(k')b(k)) = £<P' + \b\((x ~ OP + , k' ± )b x ((x + k ± )\P+) , (11) 

A 

and V is a volume factor. We are interested here in the x > £ region, since we want to 
obtain only the quark contribution to H, the only one which can be evaluated in a CQM 
with three valence, point-like quarks. 
Eq. (10) can be written: 

l r d 2 kidk + -. 

P + X 



x 5(A;+ - (x + £)P + )(tf(k + A)6(ib)) = 

= _1 /■ _£kglk^ (^__ {x + ^\ (6 t (A . + A)6(jfc)) . (12) 
2P+V 7 2 V / F^(2vr) 3 V p+ 7 

In a NR framework, states and creation and annihilation operators have to be normalized 

according to 

(P'|P> = (2tt) 3 5(P'+ - P + )5(Pi - Pi) (13) 

and 

{b(k' + , k ± ), tf(k + , k' ± )} = (2ir) 3 5(k' + - k + )5(k' ± - k ± ) (14) 

respectively. As a consequence, in order to perform a NR reduction of Eq. (12), one has to 
consider that [45]: 



|P> -> V/2P+IP) , (15) 
b(k) -> V2k+b(k + ,k ± ) , (16) 



so that in Eq. (12), in terms of the new states and fields, one has to perform the substitution 
+ A)&(fc)) = £(P' + \ b x(( x - 0P + , k ± + A±)b x ((x + k ± )\P+) - 



2y/l - eP + ^k'+2k+ J2(P'\b{(k + + A+, k ± + A ± )b x (k + , k ± ) \P) , (17) 



where use has been made of the relation 2\/l — £ 2 P + = \ / 2P + \ / 2P' + . 
Now, inserting Eq. (17) in Eq. (12), one gets 



F| + (x,e,A 2 ) = i/ 



d 2 k<dk + Jk + . \ 



X 



2 v / AFaF(2tt) 3 \P+ 
V2k'+2k+Y,(P'\ b l(k + + A+ k ± + A ± )b x (k + , k±)\P) + 0(£ 2 ) 

A 

l r d 2 k\dk + J k+ . \ 



V J (2tt) 3 \P+ 
x Y.( p '\b[(k+ + A\k ± + A ± )b x (k\ k ± )\P) + 0(e) • 



(18) 



Now, since the constituent quarks with mass m are taken to be on shell, so that ko 



+ m 2 , one has 



d 2 k ± dk^ 



(19) 



so that, from Eq. (18) and Eq. (7), one gets: 
H q (x,£,A 2 ) = JdtS^-ix + tfj x 



1 kr 



= Jdk8^t-(x + 0^ 



Y,(P'\b[(k + + A + , k ± + A ± )b x (k + , k ± ) \P) 







Ik 2 ( 


k + A) 2 Y 


I m 2 ' 


m 2 J 



+o(a= (20) 

0(f) ■ (21) 



— * — * — * 

In the last step, the definition of (non-diagonal) momentum distribution, n q (k, k + A), has 
been used, together with the fact that a NR momentum distribution describes the probability 
of finding a constituent of momentum A; in a given system up to terms of order ^ [46]. 

Summarizing, we find that, in a NR CQM, the GPD H q (x,£, A 2 ) can be calculated, for 
( 2 < 1, P < m 2 and (k + A) 2 <C m 2 (which means, in turn, A 2 <C m 2 ), through the 
following expression: 
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H q (x, £, A 2 ) = JdkS^-(x + n q (k, k + A) + U 2 , 



9 ' 9 



(22) 



The above equation, corresponding to Eq. (8) in [30], permits the calculation of H q (x, £, A 2 ) 
in any CQM, and it naturally verifies some of the properties of GPDs. In fact, the unpo- 
larized quark density, q(x), as obtained by analyzing, in IA, DIS off the nucleon (see, e.g., 
[45]), assuming that the interacting quark is on-shell, is recovered in the forward limit, where 
A 2 = £ = 0: 

q(x) = H q (x, 0,0) = J dkn q (k) 5 (x - ^ , (23) 

— # 

so that the constraint Eq. (3) is fulfilled. In the above equation, n q (k) is the momentum 
distribution of the quarks in the nucleon: 

n q (k)= J e i ^- r ' ) p q (fy)drdf f . (24) 

Besides, integrating Eq. (22) over x, one obtains 

J dxH q (x,C,,A 2 ) = J dfe^ Pq (r) , 

where p q {r) = lmv^ T ?p g (r 1 ', r) is the contribution of the quark q to the charge density. The 
r.h.s. of the above equation gives the IA definition of the charge f.f. 

J dre iK -%(r) = F"(A 2 ) , (25) 

so that, recalling that F q (A 2 ) coincides with the non relativistic limit of the Dirac f.f. 
F-?(A 2 ), Eq. (4) is verified. 

Besides, the polynomiality condition is formally fulfilled by the GPD defined in Eq. (22), 
although the present accuracy of the model, explicitly written in the latter equation, does 
not allow to really check polynomiality, due to the already mentioned effects of the Munich 
Symmetry [44]. 

The definition of H q (x,£, A 2 ) in terms of CQM wave functions can be generalized to 
other GPDs, and the relation of the latter quantities with other form factors (for example the 
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magnetic one) and other PDFs (for example the polarized quark density) can be recovered. 
Therefore the proposed scheme allows one to calculate the GPDs by using any CQM. 
With respect to Eq. (22), a few caveats are necessary. 

i) One should keep in mind that Eq. (22) is a NR result, holding for £ 2 <C 1, under the 
conditions ? « m 2 , A 2 <C m 2 . If one wants to treat more general processes, the NR limit 
should be relaxed by taking into account relativistic corrections. In this way, at the same 
time, an expression to evaluate E q (x, £, A 2 ) could be obtained. Since our main aim here is 
to describe our approach, rather than to obtain realistic estimates, we postpone to a later 
publication the discussion of a relativistic model, which will permit to study the full A 2 - 
and £-range, together with the GPD E q (x,£, A 2 ). 

ii) The Constituent Quarks are assumed to be point-like. 

iii) If use is made of a CQM, containing only constituent quarks (and also antiquarks 
in the case of mesons), only the quark (and antiquark) contribution to the GPDs can be 
evaluated, i.e., only the region x > £ (and also x < — £ for mesons) can be explored. In 
order to introduce the study of the ERBL region (— £ < x < £), so that observables like 
cross-sections, spin asymmetries and so on can be calculated, the model has to be enriched. 

iv) In actual calculations, the evaluation of Eq. (22) requires the choice of a reference 
frame. In the following, the Breit Frame will be chosen, where one has A 2 = —A 2 and, 
in the NR limit we are studying, one finds \/2P + — > M. It happens therefore that, in the 
argument of the 5 function in Eq. (22), the x variable for the valence quarks is not defined 
in its natural support, i.e. it can be larger than 1 and smaller than £. Several prescriptions 
have been proposed in the past to overcome such a difficulty in the standard PDFs case 
[33,34]. Although the support violation is small for the calculations that will be shown here, 
it has to be reported drawback of the approach. 

The issue iii) will be discussed in the next sections, by relaxing the condition ii) and 
allowing for a finite size and composite structure of the constituent quark. 
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III. GPDS IN A CONSTITUENT QUARK SCENARIO 



The procedure described in the previous section, when applied in the standard forward 
case, has been proven to be able to reproduce the gross features of PDFs [34]. In order to 
achieve a better agreement with data, the approach has to be improved. 

In a series of previous papers, it has been shown that unpolarized [37] and polarized [38] 
DIS data are consistent with a low energy scenario dominated by composite constituents 
of the nucleon. This was obtained using a simple picture of the constituent quark as a 
complex system of point-like partons, and thus constructing the forward parton distributions 
by a convolution between constituent quark momentum distributions and constituent quark 
structure functions. The latter quantities were obtained by using updated phenomenological 
information in a scenario firstly suggested by Altarelli, Cabibbo, Maiani and Petronzio 
(ACMP) already in the seventies [39]. 

Following the same idea, in this section a model for the reaction mechanism of an off- 
forward process, such as DVCS, where GPDs could be measured, will be proposed. As a 
result, a convolution formula giving the proton H q GPD in terms of a constituent quark 
off-forward momentum distribution, H qo , and of a GPD of the constituent quark q itself, 
H qoq , will be derived. 

It is assumed that the hard scattering with the virtual photon takes place on a parton of 
a spin 1/2 target, made of complex constituents. This can be the case of a spin 1/2 nucleus, 
such as 3 He, or of the proton, if this is assumed to be made of composite constituent quarks. 
The latter situation is the one we are interested in. 

The scenario we are thinking of is depicted in Fig. 1 for the special case of DVCS, in 
the handbag approximation. One parton (current quark) with momentum k, belonging to a 
given constituent of momentum p, interacts with the probe and it is afterwards reabsorbed, 
with momentum k + A, by the same constituent, without further re-scattering with the 
recoiling system of momentum Pr. We suggest here an analysis of the process which is quite 
similar to the usual IA approach to DIS off nuclei [45-47]. 



12 



In the class of frames chosen in section 2, and in addition to the kinematical variables, x 
and £, already defined, one needs a few more to describe the process. In particular, x' and 

for the "internal" target, i.e., the constituent quark, have to be introduced. The latter 
quantities can be obtained defining the "+" components of the momentum k and k + A of 
the struck parton before and after the interaction, with respect to P + and p + = \(p + p') + : 

k + = {x + i)P + = {x' + i')P + , (26) 
(k + A)+ = (x- £)P+ = (x' - £')p + . (27) 

From the above expressions, £' and x' are immediately obtained as 

A+ 

£' = ~=-r (28) 
x' = ^x (29) 

and, since ^ = — A + /(2P + ), if z = p + /P + , one also has 

^I(TtW (30) 

These expressions have been already found and used in the IA analysis of DVCS off the 
deuteron [48] and, in general, off nuclei [49]. 

In order to derive a convolution formula, a standard procedure will be adopted [45-47]. 
In Eq. (20), two complete sets of states, corresponding to the interacting constituent and 
to the recoiling system, are properly inserted to the left and right-hand sides of the quark 
operator: 

H q (x, £, A 2 ) = (P' + \ £ {\P'A)\p's')}{(P'rS' r \(p's'\} 

pi ql „/ 

E {\PrSe)\ps)H(P k S r \(ps\} \p+), 

Pr,Sr,p,s 

and since, using IA, 

{(P R S R \(ps\}\PS) = (P R S R ,ps\PS)(27r) 3 5 3 (P - P R - p)5 s ,s R s , 
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a convolution formula, valid for any GPD H q of the spin 1/2 complex target (the proton 
in the present case) in terms of the GPD H qoq of spin 1/2 structured constituents (the 
constituent quark, in the present case), go, is readily obtained: 

H q {x,£,A 2 ) = £ JdE f dpP qo (p,p + A,E)^H qoq (x',e,A 2 ) . (31) 
go J J s 

In the above equation, E is the excitation energy of the recoiling system and P qo (p, p + A, E) 

the one-body off-diagonal spectral function for the constituent quark q in the proton: 

P q MP + \^ = E(P'S\(P-p)SRAP + ^)s}((P-p)S R ,ps\PS)5(E-E R ) . 

Sr,s 

If the ^-dependence of i.e., the independence of z (cf. Eq. (30)) is disregarded in Eq. 
(31), so that the one-body off-diagonal momentum distribution 

n qo (p,p + A) = J dEP qo (p,P + A, E) (32) 

is recovered, Eq.(31) can be written in the form 

H q (x,£, A 2 ) ~ W dpn qo (p,p + A)^-H qoq (x', A 2 ) = 
go J ? 

= E jf f / dpn qo (p,p + A)5 [z - |) if,, ^, | A 2 ) . (33) 
Taking into account that 

*- | = *" [5(1 + -Z] = z + Z- + = * + £-■]£:, (34) 
Eq. (33) can also be written in the form: 

^.^^E/ 1 -^^.*. f",-,A 2 ) , (35) 

where 

tf, (z,£,A 2 ) = jdpn q MP + ^{z + ^-^j (36) 

is to be evaluated in a given CQM, according to Eq. (22), for q = uq or do, while 
H qoq (^, |, A 2 ) is the constituent quark GPD, which is still to be discussed and will be 
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modelled in the next section. One should notice that the forward limit of Eq. (35) gives 
the expression which is usually found, for the parton distribution q(x), in the IA analysis of 
unpolarized DIS off nuclei [45-47]. In fact, if in Eq. (35) the index q , labelling a constituent 
quark, is replaced by the index N, labelling a nucleon in the nucleus, in the forward limit 
one finds the well known result: 



is the light-cone momentum distribution of the nucleon N in the nucleus and gjv(^) — 
H Nq (x, 0,0) is the distribution of the quark of flavor q in the nucleon N. 

IV. A MODEL FOR THE GPDS OF THE CONSTITUENT QUARK 

The crucial problem now is the definition of H qoq (^, | , A 2 ), the constituent quark GPD, 
appearing in Eq. (35). 

As usual, we can start modelling this quantity thinking first of all to its forward limit, 
where the constituent quark parton distributions have to be recovered. As we said in the 
previous section, in a series of papers [37,38] a simple picture of the constituent quark as a 
complex system of point-like partons has been proposed, re-taking a scenario suggested by 
Altarelli, Cabibbo, Maiani and Petronzio (ACMP) [39]. 

Let us recall the main features of that idea. 

The constituent quarks are themselves composite objects whose structure functions are 
described by a set of functions <p qoq (x) that specify the number of point-like partons of type 
q which are present in the constituent of type go, with fraction x of its total momentum. 
We will hereafter call these functions, generically, the structure functions of the constituent 
quark. 




(37) 



where 




(38) 
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The functions describing the nucleon parton distributions are expressed in terms of the 
independent <p qoq (x) and of the constituent density distributions (q = u , d ) as, 

q(x, Q 2 ) = W* -*>(*, Q 2 )<p qm (-, Q 2 ) , (39) 

q J x z \ z / 

where q labels the various partons, i.e., valence quarks (u v ,d v ), sea quarks (u s ,d s , s), sea 
antiquarks (u, d, s) and gluons g. 

The different types and functional forms of the structure functions of the constituent 
quarks are derived from three very natural assumptions [39]: 

i) The point-like partons are QCD degrees of freedom, i.e. quarks, antiquarks and gluons; 

ii) Regge behavior for x — > and duality ideas; 

iii) invariance under charge conjugation and isospin. 

These considerations define in the case of the valence quarks the following structure 
function 



r(A+|) (l 



X 



,A-1 



For the sea quarks the corresponding structure function becomes, 

C 



<n = -(l-aO, (41) 



x 



and, in the case of the gluons, it is taken 

<P gg (x,Q 2 ) = -(l-xf-' . (42) 

Ob 

The last assumption of the approach relates to the scale at which the constituent quark 
structure is defined. We choose for it the so called hadronic scale /Xq [34,50]. This hypothesis 
fixes all the parameters of the approach (Eqs. (40) through (42)). The constants A, B, G and 
the ratio C/D are determined by the amount of momentum carried by the different partons, 
corresponding to a hadronic scale of /ijj = 0.34 GeV 2 , according to the parametrization of 
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[50]. C (or D) is fixed according to the value of F 2 at x = 0.01 [39], and its value is chosen 
again according to [50]. We stress that all these inputs are forced only by the updated 
phenomenology, through the 2 nd moments of PDFs. The values of the parameters obtained 
are listed in [37]. 

We, here, note that the unpolarized structure function F 2 is rather insensitive to the 
change of the sea {C , D) and gluon (B, G) parameters. 

The other ingredients appearing in Eq. (39), i.e., the density distributions for each 
constituent quark, are defined according to Eq. (22). 

Now we have to generalize this scenario to describe off-forward phenomena. Of course, 
the forward limit of our GPDs formula, Eq. (35), has to be given by Eq. (39). By taking 
the forward limit of Eq. (35), one obtains: 



so that, in order for the latter to coincide with Eq. (39), one must have H qoq (x, 0, 0) = 
<t> qoq ( x )- 

In such a way, through the ACMP prescription, the forward limit of the unknown con- 
stituent quark GPD H qoq (^, |, A 2 ) can be fixed. 

Now the off-forward behavior of the Constituent Quark GPDs has to be modelled. 

This can be done in a natural way by using the "a-Double Distributions" (DD's) language 
proposed by Radyushkin [3,43]. DD's, $(5;, a, A 2 ), are a representation of GPDs which 
automatically guarantees the polynomiality property. GPDs can be obtained from DD's after 
a proper integration. In constructing models, the DD's can represent a more appropriate 
language with respect to GPDs. In fact, the hybrid character of GPDs, which are something 
in between parton densities q(x) and distribution amplitudes 0(a), is naturally emphasized 
when the latter are obtained from DD's. The DD's do not depend on the skewedness 
parameter £; rather, they describe how the total, P, and transfer, A, momenta are shared 
between the interacting and final partons, by means of the variables x and a, respectively. 




(43) 
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As it can be argued from Fig. 2, where the DDs representation of GPDs is illustrated 
schematically, and as it is explained in [3,43], parton densities are recovered in the forward, 
A = limit, while meson distribution amplitudes are obtained in the P = limit of 
DD's. In some cases, such a transparent physical interpretation, together with the symmetry 
properties which are typical of distribution amplitudes (a — > —a symmetry), allows a direct 
modelling, already developed in [26]. 

The relation between any GPD H, defined a la Ji, for example the one we need, i.e. H qoq 
for the constituent quark target, is related to the a-DD's, which we call <& qoq (x, a, A 2 ) for 
the constituent quark, in the following way [3,43]: 

/l rl — \x\ _ 

dx / 5(x + £a- x)® qoq (x, a, A 2 )da . (44) 
-1 J-l+\x\ 

With some care, the expression above can be integrated over x and the result is explicitly 
given in [3]. The DDs fulfill the relation: 

$ qoq (x , a, A 2 ) = $ qoq (x , -a, A 2 ) , (45) 

and the polynomiality condition [2]. 

In [43], a factorized ansatz is suggested for the DD's: 

$> qoq (x, a, A 2 ) = h q (x, a, A 2 )% oq (S)F qo (A 2 ) , (46) 

with the a dependent term, h q (x, a, A 2 ), which has the character of a mesonic amplitude, 
fulfilling the relation: 

rl-\x\ 

/ h q (x,a, A 2 )da = 1 . (47) 

J-l+\x\ 

Besides, in Eq. (46) <fr qoq (x) represents the forward density and, eventually, F qo (A 2 ) the 
constituent quark form factor. 

One immediately realizes that the GPD of the constituent quark, Eq. (44), with the 
factorized form Eq. (46) and the normalization Eq. (47), fulfills the crucial constraints of 
GPDs, i.e., the forward limit, the first-moment and the polynomiality condition, the latter 
being automatically verified in the DD's description. 
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In the following, we will assume for the constituent quark GPD the above factorized 
form, so that we need to model the three functions appearing in Eq. (46), according to the 
description of the reaction mechanism we have in mind. 

For the amplitude h q , use will be made of one of the simple normalized forms suggested 
in [43], on the bases of the symmetry properties of DD's: 

c('.°)=i ( 7iy • < 48) 

Besides, since we will identify quarks for x > £/2, pairs for x < |£/2|, antiquarks for 

x < — £/2, and, since in our approach the forward densities $ qoq (x) have to be given by the 

standard $ functions of the AC MP approach, Eqs. (40)-(42), one has, for the DD of flavor 

q of the constituent quark: 

((h q (x,a)% oqv (x) + h q (x : a)^ qoqa (x))F qo (A 2 ) for x > 
% oq (x,a,A 2 ) = \ (49) 
[ -h q (-x,a)<S> qoqs (-x)F qo (A 2 ) for x < 

The above definition, due to Eq. (47), when integrated over a gives the correct limits 

[3]: 



$ qoq (x, a, A 2 = 0) \ i>0 da = <5> qoq (x) , (50) 

-l+x 

and 

fl-\x\ _ 

/ $ qoq (x,a,A 2 = 0)\x <0 da = —& qoq -(—x) . (51) 

J-l+\x\ 



Eventually, as a f.f. we will take a monopole form corresponding to a constituent quark 
size tq ~ 0.3/m: 

^o(A 2 ) = W » ( 52 ) 

1 ~-6 3 - 

a scenario strongly supported by the analysis of [40]. 

By using such a f.f. and Eq. (48), together with the standard ACMP $'s, Eqs. (40) and 
(41), in Eq. (49), and inserting the obtained <& qoq (x, a, A 2 ) into Eq. (44), the constituent 
quark GPD in the ACMP scenario can be eventually calculated. 
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V. RESULTS AND DISCUSSION 



In this section we present the results obtained for the GPD H q of the proton, for £ 2 <C 1 
and A 2 <C m 2 , according to the approach described so far. The main equation to be 
evaluated is Eq. (35), written again here below for the sake of clarity: 

qo Jx Z \Z Z J 

In the above equation, the quantity H qoq , the constituent quark GPD, is modelled ac- 
cording to the arguments described in the previous section. This means that it is obtained 
evaluating Eq. (44), where the DD of the constituent quark, & qoq (x, a, A 2 ) , is given by Eq. 
(49), calculated in turn through the f.f., Eq. (52), the function h q , Eq. (48), together with 
the standard ACMP $'s, Eqs. (40) and (41). 

The other ingredient in Eq. (35), H qo , has been evaluated according to Eq. (36): 

H qo (z, £, A 2 ) = J dpn qo (p,p + A)5^z + ^-t- 

The calculation has been performed in the Breit frame, where one has, in the NR limit 
studied, \/2P + — > M. The off-diagonal momentum distribution appearing in the formula 
above, n qo (p, p + A), defined in Eq. (32), has been evaluated within the Isgur and Karl (IK) 
model [35]. The calculation is described in [30] and the main results are listed again here 
for the reader's convenience. 

In the IK CQM [35], including contributions up to the 2hu shell, the proton state is 
given by the following admixture of states 

\N) = a s \ 2 S 1/2 )s + as'\ 2 S[/ 2 )s + a M ?Si/ 2 ) M + a- D \ i D 1/2 ) M , (53) 

where the spectroscopic notation \ 2S+1 Xj) t , with t = A,M,S being the symmetry type, 
has been used. The coefficients were determined by spectroscopic properties to be [51]: 
a s = 0.931, a s > = -0.274, a M = -0.233, a v = -0.067. 

The results for the GPD H(x, £, A 2 ), neglecting in (53) the small D-wave contribution, have 
been found to be [30]: 
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M ( 3 \ 3 ^ 2 a 2 r r 
H u (x, £, A 2 ) = 3— [-) J dk x J dk y f (k x , ky, x, £, A 2 ) 

x [/s(&x> ky, x, £, A 2 ) + f(k x , k y} x, £, A 2 ) , 



M ( 3 \ 3 / 2 a 2 r r 
H d (x, f , A 2 ) = 3— f — J J dk x J dk y f (k x , k y , x, £, A 2 ) 



x 



-^fs{k x ,k y ,x,^, A ) f(k x ,k y ,x,£,A ) 



for the flavors u and d, respectively, with 

ko 



fo(k x , k y , x, £, A 



ko + k 2 



/a(A X , k X )f a (A.y, ky) f a (A z , k Z ) , 



M 2 (x + Q 2 -(m 2 + k 2 x + k 2 y ) 
2M(x + 



fs{k x , k y , x, £, A ) — — a s + a s / 



5 _ k 2 1 A; 4 2 / A 2 - A / k 2 

6 ^ + 2^ + 3^ T + A '*J U 2 " , 



+ «M 



5 1A; 2 1A; 4 2 fc /9, 9 A9 nX r 



3a 2 \ 3 
2 



cr 



+ a s a s > 



i - 



A: 2 



2 / A 2 ^ r N 
— =- + A • k 



a 2 I 3a 2 \ 3 



f(k x ,k y ,x,£,A 2 ) = -a s a s > 



a 2 ) 3a 2 i 3 , 



— aMds' 



+ 



2P 
3a 4 



A 2 



6 ~ a 2 + 2 a* 



+ A-k 



3a 2 



A 2 -. -> 



and k = ^ m 2 + k 2 + k 2 + k 2 , being m ~ M/3 the constituent quark mass. Here we 
used the notation A; 2 = k 2 ,k = (k x , k y , k z ). 
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The harmonic oscillator parameter, a, of the IK model, can be chosen so that the ex- 
perimental r.m.s. of the proton is reproduced by the slope, at A 2 = 0, of the charge form 
factor. Such a choice is performed as follows. Integrating Eq. (35) over x, one gets: 



where F q (A 2 ) is the contribution, of the current quark of flavor q, to the proton f.f.; F qo (A 2 ) 
is the contribution, of the point-like constituent quark of flavor q , to the proton f.f.; F qoq (A 2 ) 
is the contribution, of the current quark of flavor q, to the f.f. of the composite constituent 
quark of flavor q . 

The latter is given by Eq. (52), while in the IK model one has [51]: 



with a ~ 0.015 and b ~ 0.0001. Imposing that the slope of the f.f., Eq.(61), reproduces at 
A 2 = the experimental proton r.m.s., a value of a = 1.18/ -1 is obtained. Such a form 
factor reproduces well the data at the low values of A 2 which are accessible in the present 
approach. For higher values of A 2 , it would not be realistic [51]. 

Results for the w-quark H distribution, at the scale of the model /z 2 , are shown in Figs 
3 to 5. In Fig. 3, it is shown for A 2 = —0.1 GeV 2 and £ = 0.1. One should remember that 
the present approach does not allow to estimate realistically the region —A 2 > m 2 ~ 0.1 
GeV 2 , so that we are showing here the result corresponding to the highest possible A 2 
value. Accordingly, the maximum value of the skewedness is therefore £ ~ 0.17 (cf. Eq. 
(2)), fulfilling the requirement ^ 2 < 1. The dashed curve represents what is obtained in 
the pure Isgur and Karl model, i.e., by evaluating Eq. (54). One should notice that such a 
result, obtained in a pure valence CQM, should vanish for x < £ and for x > 1. The small 
tails which are found in these forbidden regions represent the amount of support violation 
of the approach. In particular, for the shown values of A 2 and £, a violation of 2 % is found. 
In general, in the accessible region the violation is never larger than few percents. The full 
curve in Fig. 3 represents the complete result of the present approach, i.e., the evaluation 



F 9 (A 2 ) = £F^A 2 )F^(A 2 ), 



(61) 




(62) 
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of Eq. (35) following the steps and using the ingredients described in this section and in the 
previous one. A relevant contribution is found to lie in the ERBL region, in agreement with 
other estimates [4]. As already explained, the knowledge of GPDs in the ERBL region is a 
crucial prerequisite for the calculation of all the cross-sections and the observables measured 
in the processes where GPDs contribute. We notice that the ERBL region is accessed 
here, with respect to the approach of Ref. [30] which gives the dashed curve, thanks to the 
constituent structure which has been introduced by the ACMP procedure. 

In Fig. 4, special emphasis is devoted to show the ^-dependence of the results. For the 
allowed £ values, H u (x, £, A 2 ), evaluated using our main formula, Eq. (35), is shown for four 
different values of x. It is clearly seen that such a dependence is strong in the ERBL region, 
while it is rather mild in the DGLAP region, in good agreement with other estimates [4]. To 
allow for a complete view of the outcome of our approach, in Fig. 5 the x and £ dependences 
are shown together in a 3-dimensional plot. 

The results shown so far are associated with the low scale of the model, the hadronic 
scale //q, fixed to the value 0.34 GeV 2 , as discussed in Section 4. As an illustration, in 
Fig. 6 and 7 the Next to Leading Order QCD-evolution of the Non Singlet (NS), valence 
u distribution, up to a scale of Q 2 = 10 GeV 2 , is shown. One should notice that any NS 
distribution is symmetric in x, due to its definition and in agreement with the conventions 
used: 

H» s (x, £, A 2 ) = H q (x, £, A 2 ) - H q (x, £, A 2 ) = 

= H q (x,Z,A 2 ) + H q (-x,Z,A 2 ) . (63) 

In evolving the model results, the approach of Ref. [28] has been applied and a code kindly 
provided by A. Freund has been used. The evolution clearly shows a strong enhancement of 
the ERBL region. 

We have therefore developed a scheme which provides us with the GPD H in the full 
x range. This is obtained thanks to the constituent quark structure, implemented dressing 
the three quarks of a CQM, where initially only the DGLAP region of GPDs was accessible. 
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This is an important development with respect to previous work, a prerequisite for any 
attempt to calculate cross sections and asymmetries of related processes. The next step of 
our studies will be indeed to use the obtained GPDs for the evaluation of cross sections which 
have been recently measured [52,53] or will be measured soon. The estimate of cross sections 
is presently in progress, together with that of relativistic corrections, which will permit to 
enlarge the kinematical range (basically the A 2 and the related ^-ranges) where our results 
can be applied to predict or to describe the data. For the time being, the comparison with 
data of our estimates is therefore not possible. 

Anyway, as already said, the present approach fulfills several theoretical constraints. First 
of all, the forward limit, Eq. (39), provides us with a reasonable description of quark densities 
(see Ref. [37], where the IK model together with the ACMP mechanism has been applied 
to unpolarized DIS 1 ). Secondly, the x— integral of our GPD H turns out to be, formally 
and numerically, £ independent, satisfying therefore the polynomiality condition for the first 
moment, and providing us with a proton form factor, Eq. (62), in good agreement with 
the data in the low-A 2 region which is studied here. The main theoretical drawback is the 
already discussed support-violation, which in any case affects our findings, in the kinematical 
region under investigation, by a few percents at most. Another theoretical constraint which 
is satisfied by the present approach is the inequality: 



where x\ = (x + £)/(l + £) and x 2 = (x — £)/(l — £), proven in [43]. As an illustration, the 
validity of the above inequality is shown in Fig. 8 for the -u-flavor, and for A 2 = —0.1 GeV 2 
and £ = 0.1. 

Other theoretical constraints, such as the Ji sum rule [9], require the knowledge of the 
other unpolarized GPD, E, which has not yet been calculated in the present scheme. E arises 



1 In that work [37], we showed that more sophisticated quark models were producing an excellent 
description of the data [54] . 




(64) 
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naturally in a relativistic framework, where we plan to calculate it. In fact, a relativistic 
calculation will permit us to overcome the support problem and, at the same time, to enlarge 
the kinematical range of our predictions. Once a relativistic CQM is used to predict GPDs 
in the DGLAP region, the structure of the constituent quark can be introduced to access 
the ERBL region, so that the Ji sum rule and other theoretical constraints can be checked. 
Work is being carried out in that direction. 

VI. CONCLUSIONS 

Quark models have been extremely useful for understanding many features of hadrons, 
even in the DIS regime. In previous work, a thorough analysis of both polarized and unpolar- 
ized data have shown that constituent quarks cannot be considered elementary when studied 
with high-energy probes. The first feature found was the need for evolution, i.e. the con- 
stituent quarks at the hadronic scale have to be undressed by incorporating bremsstrahlung 
in order to reach the Bjorken regime, but this was not all. A second feature, we found, was 
that the constituent quarks should be endowed with soft structure in order to approach the 
data. Thus, the constituent quarks appear, when under scrutiny by high-energy probes, as 
complex systems, with a very different behavior from the current quarks of the basic theory. 
These features, which we found in structure functions, have also been recently discussed in 
form factors [40]. 

The aim of the present work has been to generalize the formalism of composite constituent 
quarks to study the generalized parton distributions. We here develop a formalism which 
expresses the hadronic GPDs in terms of constituent quarks GPDs by means of appropriate 
convolutions. In order to be able to predict experimental results we have defined a model 
which incorporates phenomenological features of various kinematical regimes. The model is 
based on Radyushkin's factorization ansatz, thus our constituent quark GPDs are defined in 
terms of the product of three functions: i) the constituent quark structure function, where 
we use the ACMP proposal [39]; ii) Radyushkin's double distributions [43]; iii) constituent 
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quark form factor as suggested in Ref. [40]. Once these GPDs are defined in this way, 
we have developed the scheme to incorporate them into any nucleon model by appropriate 
convolution. In order to show the type of predictions to which our proposal leads, we have 
used here, as an illustration, a naive model of hadron structure, namely the IK [35] model. 
However, in this latest step of our scheme, any non-relativistic (or relativized) model can be 
used to define the hadronic GPDs. 

Looking at our results, we found that the present scheme transforms a hadronic model, 
in whose original description only valence quarks appear, into one containing all kinds of 
partons (i.e., quarks, antiquarks and gluons). Moreover, the starting model produces no 
structure in the ERBL region, while after the structure of the constituent quark has been 
incorporated, it does. The completeness of the x-range, for the allowed A 2 and £, of the 
present description, is a prerequisite for the calculation of cross-sections and other observ- 
ables in a wide kinematical range. To this respect we recall that our findings hold for 
A 2 << m 2 and ( 2 « 1. Nevertheless, relativistic corrections, which permit to access a 
wider kinematical region, could be included in the approach. Our aim here has been mainly 
the illustration of our scheme, and the inclusion of relativistic corrections, together with the 
use of more sophisticated models, has been postponed and will be shown elsewhere. 

We have reminded the reader of the calculation for the diagonal structure functions and 
form factors to see how in these cases, where experimental data are available, our scheme 
leads, even with a naive quark model, to a reasonable description of the data. Thereafter, 
we have proceeded to calculate the GPDs of physical interest to guide the preparation and 
analysis of future experiments. 

This work is the continuation of an effort to construct a scheme which describes the 
properties of hadrons in different kinematical and dynamical scenarios. Our description can 
never be a substitute of Quantum Chromodynamics, but, before a solution of it can be 
found, it may serve to guide experimenters to physical processes where the theory might 
show interesting features, worthy of a more fundamental effort. 

The approach here presented can be extended to light nuclei and fragmentation functions, 
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and work is being carried out also in these directions. 
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FIGURE CAPTIONS 



Fig. 1: The handbag contribution to the DVCS process in the present approach. 
Fig. 2: Illustration of the relation between a-DDs (a) and GPDs (b). 

Fig. 3: The GPD H for the flavor u, for A 2 = —0.1 GeV 2 and £ = 0.1, at the momentum 
scale of the model. Dashed curve: result in the pure Isgur and Karl model, Eq. (54). The 
small tails which are found in the forbidden regions, x < £ and x > 1, represent the amount 
of support violation of the approach. Full curve: the complete result of the present approach, 
Eq. (35). 

Fig. 4: For the £ values which are allowed at A 2 = —0.1 GeV 2 , H u (x, £, A 2 ), evaluated using 
our main equation, Eq. (35), is shown for four different values of x, at the momentum scale 
of the model. From top to bottom, the dash-dotted line represents the GPD at x — 0.05, 
the full line at x — 0.1, the dashed line at x — 0.2, and the long-dashed line at x — 0.4. 

Fig. 5: The x and £ dependences of H u (x, £, A 2 ), for A 2 = —0.1 GeV 2 , at the momentum 
scale of the model. 

Fig. 6: The Non Singlet (NS) H GPD for valence w-quarks at £ = 0.1 and A 2 = -0.1 
GeV 2 , at the momentum scale of the model, /xjj = 0.34 GeV 2 (dashed), and after NLO-QCD 
evolution up to Q 2 = 10 GeV 2 (full). 

Fig. 7: The Non Singlet H GPD for valence w-quarks at £ = 0.1 and A 2 = -0.1 GeV 2 , 
evolved at NLO from the momentum scale of the model, /x 2 , = 0.34 GeV 2 , up to Q 2 = 10 
GeV 2 . 

Fig. 8: Full curve: the GPD H u {x, £, A 2 ), for A 2 = -0.1 GeV 2 and £ = 0.1, at the 
momentum scale of the model; dashed curve: the quantity yj 'u(xi)u(x2) / (1 — £ 2 ), where 
Xl = (x + 0/(1 + and x 2 = (x- £)/(! - f) (see text). 
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